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The light deflection of one component of a binary system due to the gravitational 
field of the other component is investigated. While this relativistic effect has not been 
observed thus far, the question arises that whether this effect becomes detectable in 
view of todays high-precision astrometry which soon will reach the microarcsecond 
level of accuracy. The effect is studied and its observability is investigated. It turns 
out, that in total there are about 10 3 binaries having orbital parameters such that 
the light deflection amounts to be at least 1 microarcsecond. Two stringent criteria 
for the orbital parameters are presented, by means of which one can easily determine 
the maximal value of light deflection effect for a given binary system. It is found, 
that for relevant binaries their orbital parameters must take rather extreme values 
in order to have a light deflection of the order of a few microarcseconds. Only in 
a very few and rather extreme binary systems the light deflection effect might be 
detectable by todays astrometry, but their existence is highly improbable. Thus, 
the detection of this subtle effect of relativity still remains a challenge for future 
astrometric missions. 
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I. INTRODUCTION 

Astrometric space missions, especiall y the ESA (European Space Agency) cornerstone 



ally 

mission Gaia, see e.g. Perryman, et al. [20011 ] . are in preparation to attain microarcsecond 



(/zas) level of accuracy in absolute positional measurements of stars and other celestial 
objects. This unprecedented accuracy of astrometric observations makes it necessary to 
account for many subtle effects which were totally negligible before. A practical model 
for as trometric observations with an accuracy of 1 //as has been formulated by Klioner 



20031 ]. where the influence of the gravitational fields inside the solar system were taken into 
account. Furthermore, a number of additional effects potentially observable at this level of 
accuracy due to various gravitational fields generated outside of the solar system were also 
briefly discussed in that investigation. One of them is the gravitational light deflection of one 
companion of a binary system in the gravitational field of the other companion (without loss 
of generality, throughout the investigation the light deflection of component B at component 
A is considered, hence, component A is considered to be the massive body, while component 
B is the light-source). This effect would change the apparent position of one component of 
a binary system. While it is clear that this effect is relatively small and even at the level of 
1 //as observable only for edge-on binary systems, the huge amount of binaries generate some 
hope that there are relevant systems where this light deflection effect becomes detectable. 
For instance, Gaia will observe 10 9 stars brighter than 20 th apparent magnitude. Detailed 
numerical simulations predict the detection of about 10 s (resolved, astr ometr ic, eclipsing, 



spectroscopic) binary systems by Gaia mission, see Zwitter & Munari (2004| . which is a 



considerable increase comp ared to the 10 binary systems known so far, see " Washington 



Double Star Catalog 1 ' [2011 



Let us consider this argument a bit more quantitatively. For a simple estimate of the 
expected order of magnitude in lig ht deflection, the cla ssical lens equation (in the for m give n 
by Eq.(67) in Fritelli, et al. |200o| . Eq. (24) in Bozza |2008| | or Eq. (23) in Zschocke [201 1| ) 



is applied. In terms of orbital elements of a binary system it can be written as follows: 




■ , A 2 „ m A A . , 
V 9 = 7; \ \ cos 1 + 16 sin i cos z . (1) 

Here, ip is the light deflectin angle, % is the inclination, A is the semi-major axis and r is the 
distance of center-of-mass of the binary system from the observer, and the Schwarzschild 

radius is m = — — , where G is the gravitational constant and c is the speed of light, and M 

c 2 

is the stellar mass of component A of the binary system. From (TjQ) one obtains the maximal 
possible light deflection for edge-on binaries, i.e. attained for the case when the inclination 
is exactly 90°: 



/ M A pc 

Here, AU = 1.496 x 10 11 m is the astronomical unit, pc = 3.086 x 10 16 m stands for parallax of 
one arcsecond, and M & is the solar mass. According to this formula, the choice of moderate 
values like M = M Q and A ~ 100 AU would result into significant light deflection effect 
on microarcsecond level even at large distances of about r ~ 100 pc. A meaningful value 
for the density of stars in the solar neighborhood, 0.025 binaries pc -3 , implies already 10 5 
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binaries inside a sphere of r = 100 pc. Thus, one might conclude among them there are a 
few relevant edge-on binary systems. 

However, in order to estimate quantitatively the number of relevant systems, one needs 
to know the probability for the occurrence of such edge-on binaries which have a given light 
deflection depending on their orbital parameters like inclination, mass, distance and semi- 
major axis. Such a relation between a given light deflection and orbital parameters is given 
by a so-called inclination formula. Since the inclinations of binary systems are of course 
randomly distributed, it is meaningful to resolve su ch an inclination formula in terms of 
inclination. As it has been shown by Klioner, et al. j2003 ] (see Appendix |D] for some basic 



steps), such an inclination formula can be obtained by means of the analytical solution of 
light deflection in standard post-Newtonian approach, and is given as follows: 



7T 

i 

2 



/ _ M /ias pc\ . , 

< 2 arctan 0.0197 — — — . 3 

kms V M Q ip r J 



For a better illustration of the inclination formula, relation is rewritten in terms of 
angular degrees instead of radians: 

o M /ias pc , . 

90° -i < 2.25° — — — , 4 

where also arctan x = x + 0(x 3 ) has been used. According to this relation, the inclination 
i of a binary system with stellar mass M and at distance r must not deviate from the 
edge-on value 90° too much in order to have a given light deflection ip. For example, for 
a hypothetical binary star with M = M Q situated at a distance of r = 10 pc the light 
deflection effect attains 1 /ias only if | 90° — % \ < 0.225°, which means that the probability 
to observe this binary at a favorable inclination is only about 0.2%. 

But even this estimation is still much too optimistic. As a concrete example of todays 
high-precision astrometry, let us consider one important parameter about the astrometric 
accuracy of the Gaia mission: the accuracy of one individual positional measurement in the 
most ideal case (bright star, i.e. 10 th magnitude) amounts to be 25 /ias, which implies ip > 
25 /ias. Furthermore, inside a sphere of 10 pc around the Sun almost every star and binary 
system is known already by the data of Research Consortium on Nearby Stars (RECONS) 
2012j |. Since inside that sphere there is no binary system having a light deflection on 



microarcsecond level, one has to take at least r > 10 pc. By taking into account these both 

remarks, one obtains = in relation (HI). Therefore, even in the best case one has 

cp r 250 

to conclude | 90° — i \ < 0.01° M/M Q , that means (for M = M & ), the probability to observe 
such a binary at a favorable inclination is practically only about 0.01%. And the binaries 
must be, in fact, almost edge-on in order to have a light deflection which can be detected by 
todays astrometry. Accordingly, while relation ([2]) triggers the hope about the existence of 
many relevant binary systems, from relation (J3J) one concludes that the number of relevant 
binary systems is considerably reduced. 

An estimation of the total amount of binaries depends on many different parameters, 
like mass, semi-major axis, inclination and distances of the binaries. Therefore, a simple 
estimation is not so straightforward as one might believe. Moreover, the relation (T5]) has 
been obtained with the aid of classical lens equation, while relation fll]) has been obtain by 
means of standard post-Newtonian approach. These both approaches have different regions 
of validity. However, a rigorous treatment of the problem of light deflection in binary systems 
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implies the need of an analytical formula which is valid for such kind of extreme astrometric 
configurations like t he bin ary systems are. Recently, a generalized lens equation has been 
derived by Zschocke |201ll | , which allows to determine the light deflection of binary systems 
on microarcsecond level. One aim of this study is, therefore, to reobtain the criteria fl2]) and 
fll]) as stringent conditions from one and the same approach, i.e. with the aid of generalized 
lens equation. This is possible, because in the corresponding limits the generalized lens 
equation agrees with the classical lens equation and the standard post-Newtonian solution. 
Another aim is, to derive an inclination formula like fl3]) for binary systems which allows to 
determine the needed inclination for a given light deflection angle and as a function of the 
orbital parameters. For that one has to take into account the distribution of stellar masses 
and the distribution of semi-major axes in binary systems. Finally, the aim of this study is, to 
determine the total number of relevant binaries having a light deflection on microarcsecond 
level, and to investigate the possibility to detect this effect of light deflection by todays high 
precision astrometry. 

The article is organized as follows: In Section [TT] some basics about orbital elements of 
binary systems are given. The generalized lens equation and the inclination formula are 
presented in Section IHIl In Section IIVI two stringent conditions on the orbital elements of 
binary systems (astrometric, spectroscopic, eclipsing and resolved binaries) are presented, 
which allow to determine whether or not the binary system will have a light deflection of 
a given magnitude. The total number of binaries which have a given light deflection for 
an infinite time of observation is estimated in Section |V] while the more practical case of a 
finite time of observation is considered in Section I VI I The special case of resolved binaries is 
considered in Section IVHI For that, the specific instrumentation of Gaia mission is considered 
in some detail as the most modern astrometric mission with todays highest possible accuracy. 
A Summary is given in Section IVIIII 



II. ORBITAL ELEMENTS OF A BINARY SYSTEM 

Consider a binary system, component A with mass Ma at coordinate ta and component 
B with mass Mb at coordinate r B - In order to express the light deflection effect in terms 
of orbital elements, spherical coordinates are introduced, illustrated by FIG. [TJ 

The center of coordinate system is located at the CMS, i.e. 



^CMS 



(M A v a + M B r B ) 



Ma + Mb 

Thus, the vector r, which points from CMS to the observer, is given by 

/ r cos u sin i \ 

r sin a; sini 



(5) 



V 



(6) 



r cos i 



J 



where r = \r\, the argument of periapsis is denoted by u, and % is the inclination, see 
FIG. [TJ The solution of equation of motion yields for vectors va and vb the expression 
given by Eqs. (IA22I) - flA26j) . The vector X\ points from the mass center of massive body to 
the observer, and vector Xq points from the mass center of massive body to the source, see 
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observer 




FIG. 1: Seven orbital elements which define the orbit of a binary system: distance vector r, 
semi-major axis A, inclination i, eccentricity e, eccentric anomaly E, periapsis u, and mass ration 
Ma/Mb- The orbit of the binary system spans the (x,y)-plane and the z-axis is perpendicular 
to the orbital plane. The x-axis is oriented along the semi-major axis of the orbit of the binary 
system, while the y-axis is perpendicular to the x-axis. The vector r is directed from the center-of- 
mass (CMS) of the binary system, see Eq. ([5]), to the observer. The center of spherical coordinate 
system is located at the CMS of binary system, i.e. t*cms = 0- The inclination < i < tt is the 
angle between r and z-axis; i = ir/2 is called edge-on and i > ir/2 corresponds to retrograde orbit. 
The dotted line indicates the projection of r onto orbital (x, y)-plane, i.e. z-component of r equals 
zero. The angle between this projection and x-axis is called argument of periapsis < u < it. The 
orbital elements semi-major axis A, eccentricity < e < 1 and mass ratio Ma/ Mb govern uniquely 
the geometric shape of both ellipses. The eccentric anomaly < E < 2tt (not plotted here), is 
defined in Eq. ()A16|) of Appendix [X] and determines the actual position of the bodies A and B on 
their orbit. 

also FIG. [TJ The coordinates of these vectors can be expressed by the orbital elements of 
the binary star as follows: 



/ 




r cos u sin i — 



Ma 



1 + 



M B 



X\ = r — va = 



Ay/l-e 2 smE 



(7) 



r sin u sin i — 




\ 



r cos i 



I 
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/ A {cosE -e) \ 



x = r B - r A 



A vi - e 2 sin £ 



(8) 



V o / 

Here, A is the semi-major axis, e is the eccentricity, and E is the eccentric anomaly, see 
Appendix |A] The vectors ([7j) and will be used to express the light deflection in terms of 
orbital elements of the binary system. 



III. INCLINATION FORMULA FROM GENERALIZED LENS EQUATION 

A scheme of light propagation of a signal emitted at component B in the gravitational 
field of component A is shown in FIG. [2J The vector X\ points from the mass center of 
massive body to the observer, and vector x Q points from the mass center of massive body 
to the source, and we define R = xi — x , the absolute value R — \R\ and unit vector by 
k = R/R. Furthermore, the impact vector d = k x {x\ x k) is defined, and its absolute 
value is denoted by d = \d\. The Schwarzschild radius of massive body, i.e. of component 

A of the binary system is denoted by m = — — . 

c 2 



light path 




component A (deflecting body) 



FIG. 2: Binary star composed of component A being the massive body, and component B consid- 
ered to be the light-source. 

For determining the light deflection in weak gravitational field s the re are two essential 
approaches: standard post-Newtonian ap proac h, e.g. Brumberg )l99l| . and classical lens 
equation e.g. Eq.(67) in Fritelli, et al. |2000l ]. Eq. (24) in Bozza |200«j | . or Eq. (23) in 



Zschocke [20111 ] . While the first approach is restricted by the condition m 3> d, the second 
approach is only valid for the case that source and observer are far from the massive body, 
especially for a = k ■ x\ ^> d and b = —k ■ xq ^> d (for geometrical illustration of a, b and 
d see FIG. [2]). However, in binary systems extreme configurations are possible like d = or 
6 = 0. Therefore, in order to investigate the light deflection in binary systems one needs a 
generalized lens equation which is valid in such extreme configurations where the standard 
post-New tonian approach as well as the classical lens equation cannot be applied. Recently, 



Zschocke has derived a generalized lens equation, which allows to determine the light 



deflection in such extreme astrometric configurations like the binary systems are: 

1 / [dP m Xq X] — Xq ■ X\ d\ 
f = oh/-2+8- D _° ~- ■ (9) 
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Actually, the lens equation has two solutions, but here only one solution is considered, while 
the second solution represents just the second image of one and the same source which is not 
relevant in our investigation. The generalized lens equation is valid up to terms of the order 

and the absolute value of their total sum can be shown to be smaller or equal to 




15 7T m 2 „ L 



— Here, d' = — is Chandrasekhar's impact parameter, see Chandrasekhar [1983 



4 d' 1 1 E 

where L being the orbital momentum and E is the energy of the photon in the gravitational 

field of massive body. Basically, the light-ray of component B cannot be observed if d' 

is smaller than the radius of massive body A. For stars, the radius is much larger than 

2 

Tft 

Schwarzschild radius m, hence <C 1. Furthermore, the generalized lens equation (JH]) is 

finite for d — > and b = —k-x — > 0, both of whi ch are possible astrometric configurations in 
binary systems. Furthermore, in Zschocke 201l| it has been shown that the generalized lens 
equation fl9]) yields in the appropriate limits the correct standard post-Newtonian solution 
and the classical lens equation, hence provides a bridge between these essential approaches. 



2.0 



^ 1.0 




0.5 - 



0.00 1 - — 1 J — — 

0.0 1/2 71 7C 3/2 7C 2 7C 

E [rad] 

FIG. 3: Typical light curve of a binary system, determined using generalized lens equation or 
(IBlip . respectively. The parameters chosen are: distance r = lpc, semi-major axis A = 100 AU 
inclination i = jj^ ir, mass Ma = 2Mq, mass ratio = 2.0, eccentricity e = 0.25, argument of 
periapsis w = ^. 

In the following, Eq. fl9]) is applied in order to determine the light deflection in binary 
systems. For that, the coordinates x and X\ are used in the form as given by Eqs. (J7|) and 
([H]), respectively. A typical light-curve of a binary system, calculated by means of generalized 
lens equation is shown in FIG. [21 

Now, an inclination formula is derived from the generalized lens equation Q. The impact 
of eccentricity is neglected, circular orbits e = are considered, implying u = 0. Thus, the 
coordinates Xq and X\ are simplified to the expressions fIBlj) and (IB2I) given in Appendix IB1 
Furthermore, the maximal value of light deflection of a binary system is of interest, i.e. the 
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astrometric configuration E = is considered here. Then, by inserting these coordinates in 
the generalized lens equation one obtains (see Eq. (1B14[) in Appendix |B]) up to terms of 

A / Tn A 
the order O I — \ 

T \ T T 



M 2 . m A . , A , 
— cos 2 1 + 8 (1 + sin i) | cos i 

£t iy ry ry 



(10) 



The minimal value <£> m i n = <p (i = 0) and maximal value Wax = V 9 (* = f ) °f light deflection 
for the astrometric position E = follow from Eq. ( TTUj) : 



V-'mir 




A 2 



771 



— + 8 ^2 — 



M pc 
0.0197 uas — — 

M r 



\JmA 



M A pc 
200/^as W — — — 
P VM AU r 



(11) 
(12) 



777 



where in (111!) terms of the order (9 — - have been neglected. The expression ffTUl) can be 

yr Ay 

reconverted in terms of inclination (see Appendix IC]) : 



7T 

2 



arccos 



P 



+ 



(13) 



where 



P 



] m 2 A — 4 m r 2 ip 2 
A (r 2 y? 2 + 4 m 2 ) 

A 2 r 2 (y9 2 + 4 777 A r 2 Lp 2 - 4 7T7 2 A 2 



4 4 



A 2 (r 2 ip 2 + Am 2 ) 



(14) 
(15) 

The inclination formula f|T3|) yields the upper limit for i of a binary system in order 

to have a given value of light deflection (p. Note, that the values of ip cannot be chosen 
arbitrarily, but they are restricted by ip m i n and y9 max given by Eqs. (fTTl) and f|T2|) . respectively. 

The inclination formula f|T3|) can considerably be simplified. From (|14p and f lT5|) one 
obtains by series expansion 

P = - 4 - + 8 ZFT? + (™ ) , (16) 



r 2 (p 2 



m 
4-r + 



m 



A r 2 ip 2 



i 11 
A m r ip _ / ,\ 

4 ^ + ^ +0 M- 



(17) 



7» 



By means of f|T2|) . the last term in (|17jl can be estimated to be smaller than 4 — <C 1. Here, 

it should be underlined that ^ <C -P^j even at large distances r ~ 10 3 pc and small values 
for semi-major axis A ~ 1 AU. Thus, one obtains 



arccos 1 — 



777 

r 2 </) 2 



2 arctan 2 



77 i 



(18) 
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(7T1 3 \ / 772 \ 

— — - and Of — ); here the relation arccosfl — 8x 2 ^ 
r 3 (f 3 J \AJ 

2 arctan2x + (9 (x 3 ) for i<1 has been used. It should be underlined, that the applicability 
of ( I18p is restricted by the condition (1D5|) given in Appendix |D] and by d ^> m. Here, it 
should be noticed that x = 2 -P^j <C 1, even in such an extreme case like r = 1 pc, m = m Q 
and (p ~ 1 /ias one obtains a small number x = 0.019; for an analytical proof use the exact 
expression for y? m in- Due to ^ jr^, the impact of semi-major axis is of lower order and 
can be neglected in the inclination formula. The inclination formula f lTB]) can be expressed 
in terms of dimensionless quantities as follows: 



CM uas pc\ , . 

0.0197 — — — 19 
M <p r 1 ; 



Note, that express ion (] T9l) agrees with an inclination formula ([3]) derived at the first time 
by Klioner, et al. |2003| : the arguments of their work are represented in Appendix [D] The 
simplified inclination formula f fl9|) is not only useful for straightforward estimations about 
the order of magnitude, but it's simple structure provides also an obvious comprehension 
about the interplay of the individual terms. 



IV. STRINGENT CONDITIONS ON ORBITAL PARAMETERS FOR BINARY 

SYSTEMS 

In this Section, two stringent conditions on the orbital elements are highlighted for bina- 
ries having a given light deflection. These strict conditions are valid for any binary system: 
astrometric, spectroscopic, eclipsing and resolved binaries. 

The first stringent condition follows from the maximal light deflection angle (Tl2|) . given 

by 

[M A~ pc . . 

^""WmJauT' (20) 

It represents a strict criterion for the maximal light deflection of a binary system with given 
Schwarzschild radius of component A, given semi-major axis A, and given distance r between 
binary system and the observer. 

The second stringent condition on the orbital elements follows from the inclination for- 
mula in the simplified form as given by Eq. ([19]) . For a better illustration, this condition is 
given in terms of angular degrees instead of radians. Using arctanx = x + 0(x 3 ) one obtains 

„ M uas pc 

90° -i \< 2.25° — — — . 21 

According to this strict condition, the inclination i of a binary system with mass M (recall 
M is the stellar mass of component A) and at distance r must not exceed the given value 
in order to have a light deflection ip. 

Both these stringent conditions Eqs. (120]) and (12T]) were already stated in the introductory 
Section by Eqs. (|2]) and (j4j), respectively. However, is should be underlined here, that 
both ( 120]) and (I2ip were obtained with the aid of one and the same approach, namely the 
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generalized lens equation, while d2J) and (jlj) were obtained by means of the classical lens 
equation and the post-Newtonian solution, that means by two different approaches. 

The observability of light deflection effect in binaries implies the realization of both these 
stringent conditions (1201 and fj2~T]) simultaneously for a given binary system. But even if 
a given binary system satisfies these both conditions, the observability of light deflection 
effect is not guaranteed, because the astrometric position E = has to be reached during 
the time of observation. Nonetheless, as soon as the orbital elements r, A, m and % of the 
binaries are known, these both stringent conditions (120]) and ( 12 ip allow to find a possible 
candidate for being a relevant binary system for a given light deflection (p depending on the 
instrumentation of the observer. However, as it will be shown in Section |V] and Section IVIj 
the existence of such systems is highly improbable. 



V. TOTAL NUMBER OF BINARIES WITH A GIVEN LIGHT DEFLECTION 
FOR INFINITE TIME OF OBSERVATION 



In the previous Section, the conditions on orbital parameters for a binary system have 
been determined in order to have a given magnitude of light deflection (p. In this Section, 
the total number of such relevant binaries is determined. In order to estimate the total 
number of binaries having a given light deflection (p, the following formula is applied: 



Rtx 



N(<p) = J d 3 rp(r) J dAf(A) J dp /(/*) P(i) . (22) 

-Rmin j4min A'min 

Here, p(r) is the density of binaries, f(A) is the semi-major axis distribution of binary 

M 

systems, /(//) is the mass distribution of stars where p = — — is the mass-ratio of massive 

body (component A) and solar mass. The probability distribution P(i) to find a binary 
system with given inclination < % < it, is a function of distance r, semi-major axis A, mass 
ratio /i, and the given light deflection angle ip. According to (TT3"]) . the probability distribution 
P(i) is given by (the inclination of binary systems is of course a random distribution) 

P(i) = I arccos ( -| + , (23) 

where p and q are given by Eqs. (fl4"|) and (|T5l) . 

For the minimal distance of a binary system from the Sun one can safely assume -R m i n = 
1 pc. From (fl2l) follows that beyond a sphere with radius -R max = 2000 pc only a very few 
exceptional binary systems might have a light deflection of more than 1 //as. The Sun is 
located in the Orion arm which has about 1000 pc across and approximately 3000 pc in 
length. For the estimate according to (122]) . it is meaningful to assume that the stars are 
homogeneously distributed inside the Orion arm. For the unifor m sta r density p s tars = 
0.1 star pc~ 3 , a value which is in line with the data of RECONS [2012j |. Furthermore, a 



common presumption is, that about 50 per cent of all stars are compone nts of a binary or 



multiple system, see Duquennoy & Mayor [19911 ] and Halbwachs, et al. |2003j . Then one 
obtains for the density of binaries 

p(r) — 0.025 binaries pc~ 3 . (24) 
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Let us now consider the distribution of semi-major axis A in binary systems. Statistical 
investigations show that the distribution of binary semi-major axis is flat in a logarithmic 
scale over the range of six orders of magnitude, that means assumed to be valid in the large 
range A min = 10 R Q < A < 10 4 AU = A max , see Kouwenhoven & de Grijs [2008]. The lower 



limit A m i n is determined by the semi-major axis at which Roche lobe overflow occurs, while 
the upper limit A max depends on how large the averaged star density is. The logarithmic 
distribution is known as " Opik's lavJ 1 (1924) named after its discoverer, and is given by 

f(A) ~ — , a law which has also been confirmed by recent investigations, see Poveda, et 



A 



al. [20061 ] . This distribution is a consequence of the process of star formation as well as of 
the dynamical history of the binary system, and one can take this law as a given fact for 
numerical studies. Accordingly, see Appendix [Ej 

f(A) = - A (m^)" 1 - (25) 

min / 



A V A 



Furthermore, for the mass distribution / (/i) let us recall the initial mass function (IMF) 
which is the probablity that a star is newly formed with a stellar mass M and is frequently 
assumed to b e a power law f(M) ~ M~ a . Originally, the IMF has been introduced by 



Salpeter [19551 ] for solar neighborhood region, who assumed the value a = 2.35 and a validity 
region for stars with masses between 0.4 M & and 10 M & . During the past decades the 
IMF has been refined by several investigations. Especially, the numerical values of slope 
param eter a and regions of valid ity have been proposed in subs equent inve stigat ions, e.g. 
Scalo 1986!]. Robin, et al. [2003 ] . Nin kovic & Trajkovska [20061 ] . Ninkovic |l995 |. Kroupa 



2001]; for a review see Kroupa [2002j ]. Moreover, the IMF does not necessarily coincide 
with the real mass distribution of stars, because IMF describes mass distribution of a star 
formation, while the solar neighborhood mainly consists of evolved stars of main sequence. 
Here, for simplicity this distribution is used as a given fact with a = 2.35, and the proposed 
region of validity /i min = 0.4, and /i max = 10 is assumed. According to IMF, one finds for 
Q^l (see Appendix [E]) 

/imax /i m in 

In ord er to motivate that distribution further, one can compare fT26|) with the RECONS data 
2012j ] where one finds a fair agreement. Using fT23|) - fT26|) . the results of the estimate fT22|) 



are shown in FIG. HI recall, when evaluating this integral one has to take into account the 
boundary conditions given by Eqs. ( TTTj) and (fT2"|) . 

According to FIG. HI in total there are about iV ~ 10 3 binaries having a light deflection 
of at least <p = 1 fias. Here, the observation time is assumed to be infinite, so that it is 
guaranteed that the optimal configuration E = is reached. 



VI. TOTAL NUMBER OF BINARIES WITH A GIVEN LIGHT DEFLECTION 

FOR FINITE TIME OF OBSERVATION 

In Eq. f f22|) the number of binaries with a given maximal possible light deflection (p has 
been determined, just by taking for eccentric anomaly the value E — 0, that means the 
ideal configuration where the light deflection takes its maximal value (note, the eccentricity 
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FIG. 4: Total number of binaries according to Eq. (|22|) . having parameters such that the light 
deflection of component B at component A is larger than a given value for cp. Note, that for one 
individual observation the astrometric accuracy of Gaia is about 25 pas and the end-of-mission 
accuracy is about 5 pas in the most ideal case (bright star). The observation time is assumed to 
be infinite, so that it is guaranteed that the optimal configuration E = is reached. 



e = 0). It is, however, obvious that during the most part of the orbital motion one will 
have E ^ and the light deflection will be much smaller than the maximal possible light 
deflection angle (p. On the other side, the orbital period T of relevant binaries, given by 
Eq. (IA2ip . will easily exceed the time of observation; for instance, the Gaia mission time is 
about T mission ~ 5 years. Therefore, it will be not very probable, that the component B will 
be just at the relevant position near the value E = 0, where the light deflection becomes 
observable on microarcsecond level. In order to determine that number of observable relevant 
binaries, one has to extend Eq. (I22p as follows, 



N(<p) 



J d 3 r p(r) 

^min ^min 



pmax 

dAf(A) J dfif(ii)P(i)P(E). 



(27) 



Here, P(E) is the probability for the binary system to be in the region E, where the light 
deflection is larger than a given value for (p. 

In the very same way, as applied for the derivation of the inclination formula (TIB"]) , one 
can reconvert (IB 15}) in terms of eccentric anomaly E and one finds the eccentric anomaly 
formula: 



E 



± arccos 




(28) 



where p and q are given by Eqs. (114]) and (115]) . For a given value of light deflection (p, the 
formula (12 8p yields the value of eccentric anomaly E of a binary system characterized by 
semi-major axis A and Schwarzschild radius (or mass) m at a distance r. However, the 
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FIG. 5: Total number of binaries according to Eq. (|27|) where the binary system reaches the optimal 
configuration E = during an assumed observation time of T m ; ss i on = 5 years, and having orbital 
parameters such that the light deflection of component B at component A is larger than a given 
value for ip. Note, that for one individual observation the astrometric accuracy of Gaia is about 
25 yuas and the end-of-mission accuracy is about 5/ias in the most ideal case (bright star). 



values of (p cannot be chosen arbitrarily, instead they are restricted by y? r 
and v?ma 



<p[E 



71 
±- 

2 

71 



= (p (E = 0) given by (of course, only astrometric positions with < E < — are 

71 

taken into account, because for the area — < E < n the light deflection is negligible): 



Pi 




A 2 



m A A 



m 



— + 8 ^2- 



M pc 
0.0197 uas — — 

P M r 



ymA 



M A pc 
200/xas J— — — 
p VM Q AU r 



(29) 
(30) 



where in 



m 



terms of the order O — - have been neglected. These expressions resemble 

V r A I 



the corresponding expressions in Eqs. ( ITT]) and (fT2l) . According to Eq. ( |28|) . the region where 
the binary system has a light deflection larger or equal (p is given by 2 E. One has also to 
take into account, that during Gaia mission time T mission the component B moves along the 
orbit and could move into the region IE. Therefore, the probability P(E) that the binary 
system is during Gaia mission time at least ones inside the relevant astrometric position 
with the value E in (I2"gj) . is given by 



P{E) = Vi 



1 



arccos 



71 




P \ , -^mission 



(31) 
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where the operator V\ is defined by 

! f f (32) 

The probability distribution (|3T|) has to be implemented in Eq. (127)1 in order to determine the 
number of binary systems having a given light deflection (p and to be observable during Gaia 
mission time T mission . The results of Eq. (|27j) are shown in FIG. [5j recall, when evaluating 
this integral one has to take into account the boundary conditions given by Eqs. (fTTj) and 
(TT2"|) . Accordingly, there are only a very few binaries ~ 10 2 having a light deflection of at 
least ip = 1 //as, if the astrometric position E = is reached during an assumed observation 
time of T mission = 5 years. 



VII. SPECIAL CASE: CONDITIONS ON ORBITAL PARAMETERS FOR 

RESOLVED BINARIES 

In this Section, the special case of a resolved binary system is considered. In order to 
investigate the detectability of the light deflection effect in binaries, todays most modern 
astrometric mission, the ESA cornerstone mission Gaia, and its instrumentation which pro- 
vides the highest possible astrometric accuracy at the moment, will be considered here as a 
concrete example. 



A. Resolving power of Gaia 



The core of Gaia optical instrumentation consists of two identical mirror telescopes, 
ASTRO-1 and ASTRO-2, with a rectangular pupil whose dimensions are A = 0.50 m, B = 
1.4 5 m, a nd / = 35 m is the eff ective focal length. The intensity is given by, see Lattanzi, et 
al. [l998| and Lindgren [l998| : 



/<»*,*) = /.(^^), (33) 

where za = n A/AsuiOa, z-q = h B/AsuiOb, A and B the width and length of rectangular 
mirror, A is the wavelength of incident light-ray, and 0a, 9b are the angle of observation, i.e. 
the angle between the axis of the rectangular aperture and the line between aperture center 
and observation point. The intensity of incident light-ray at 9a = 0, 9b = is denoted 
by I . The function I(za,zb)/Iq in Eq. ( 133]) is the (by 9a = 9b = normalized) Point 
Spread Function (PSF) for monochromatic incident light with wavelength A for a rectangular 
aperture. The optical spectrum of stars is A = (350 — 750) nm. In FIG. E]the PSF for an 
incident monochromatic light-ray with A = 350 nm is represented for Gaia telescopes. 
Most of the light is concentrated in the central bright rectangular shaped pattern. The 
length /a and width Zb of this rectangle is determined by the first zero-roots of (1331) at 
za — 7T and zb — 7T, respectively. From that it follows, that sui9a = vr X/(n A) = A/A, and 
sin9B = vr A/(7r B) = A/B. Furthermore, if the diffraction pattern is sho wn on a screen at 
a distance /, then the length and width is given by, see Hog, et al. [1997[: 



^a = 2^, (34) 



16 




FIG. 6: Point Spread Function (PSF) for a rectangular telescope according to Eq. ([33]). The inci- 
dent monochromatic light-ray has a wavelength of A = 350 nm. The parameters of the rectangular 
telescope are: A = 0.5 m, B = 1.45 m. 



£b = 2^, (35) 

where / is the focal length of the optic, i.e. of the rectangular Gaia mirror. Accordingly, the 
given numerical values A = 0.50 m, B = 1.45 m, and A = 350 nm result in La = 49.0 /im and 
Lb = 16.9 fim for the length and width of "Airy rectangle" of Gaia optics. Note, that the 
"Airy rectangle" has the same order of magnitude than the pixel size (10/im x 30yum) of the 
110 CCDs (Charge-Coupled Device) sensors of astrometric field part of the focal plane. In 
order to separate two pointlike sources, the distance between their centers of the rectangle 
has to be larger than either La or Lb- Since Lb < La, in this study the better resolution 
value L B is used, which corresponds to a resolution angle of 

' = £ = (36) 

The resolving power is the minimal angular distance between two objects to get separable 
by Gaia instrumentation. With the parameters given above one obtains the resolving power 
5 of Gaia optics: 

5 = 0.24 x 10" 6 rad = 49.7 mas . (37) 
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In what follows this parameter is of fundamental importance in order to determine the ability 
of Gaia to determine the light deflection in binary systems. 



B. Orbital Parameters of Resolved Binaries Observable by Gaia 

In this Section, the question is addressed which and how many binary systems can be 
separated by Gaia instrumentation among all those relevant binaries found in the previ- 
ous Section; see FIG. El In average, Gaia will observe each object 80 times, but will not 
constantly observe these objects during mission time. However, for simplicity the scanning 
law of Gaia is approximated by assuming a permanent observation of all objects during the 
whole mission time. 

Furthermore, visual binaries are considered in this Section, i.e. binaries which are sep- 
arable by Gaia telescopes. The both largest telescopes of Gaia have a resolution angle S 
discussed in the previous Section, see Eqs. fl36l) and fl37]) . For binary systems, this resolution 
angle 5 corresponds to a minimal distance between the components A and B to get separable 
within Gaia optics. Using ([B6]) and f )36|) one obtains the condition 

d = A |cosz| > 5r = — r, (38) 

B 

where r is the distance between the remote objects and Gaia observer. This condition is by 
far much more important than taking into account the effect of finite radius of the stars, 
which would imply A |cosi| > Ra, where Ra being the radius of component A. By inserting 
the extreme case A |cosi| = 5r into ffTU]) . one obtains 




.... , . m A \ m A 1 . . 

<P = o h/^ + 8(l + sin*) ---^1 (39) 

where sinz ~ 1 has been used, and terms of higher order O (m 2 ) are neglected. Relation 
( |39|) is an expression for the maximal light deflection angle of a binary system when taking 
into account the resolving power of Gaia. Eq. (139]) is a much stricter restriction than the 
generalized lens equation f TTO]) . because fl39~l) determines the light deflection angle only of 
those binary systems having a resolution angle 5 of Gaia optical instrumentation, while ( TTTTj) 
determines the light deflection angle of any possible binary system. Note, that from Eq. (13"9"1) 
follows the maximal possible distance of visual binaries for a given light deflection: 




where in the last expression the optimal values 5 = 0.24 x 10~ 6 rad and ip = 25 /ias has been 
used. The condition (j40j) can also be written by 



A>30A UT f — . (41) 



M Q r 2 

~W p^ 2 

These both conditions fj4*0"l) and ([4TJ) imply rather extreme orbital parameters on visual 
binaries. For instance, condition fj4*0~l) implies a maximal distance of r < 18 pc for solar-mass 
type binaries even with a huge semi-major axis of A = 10 4 AU, while condition (T4T]) implies 
a large semi-major axis for solar-mass type binaries outside a sphere of r > 10 pc. It is 
almost certain, that such extreme parameters will not be realized in nature. 
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VIII. SUMMARY 



In this study, the light deflection in binary systems has been considered. While there 
is absolutely not any doubt about the existence of this relativistic effect, is has not been 
observed so far. To investigate this effect of light deflection, an inclination formula flT3|) ha s 



been derived by means of generalized lens equation ([9]) obtained recently by Zschocke [2011 



and these both equations are the theoretical basis for investigating the light deflection effect 
in binary systems. A simplified inclination formula has been presented by Eq. (Tl9|) and its 
validity has been discussed in som e detail. This simplified inclination formula has also been 
obtained by Klioner, et al. [20031 ] by an independent approach. Furthermore, two stringent 
conditions on the orbital parameters have been given by Eqs. (1201) and (1211) . These both 
stringent conditions are relations between the orbital elements of a (resolved, astrometric, 
eclipsing, spectroscopic) binary system for a given magnitude of light deflection, and allow 
to find a relevant binary system in a straightforward way. 

In Section |V] the total number of binaries with a given light deflection has been deter- 
mined by means of the semi- major axis distribution according to " Opik's law" and the mass 
distribution according to " Salpeter's mass distribution" . Since the inclinations are randomly 
distributed, the inclination formula allows to estimate the total number of relevant binaries 
with the aid of Eq. f[2"2"j) . It turns out, that in total there exist about 10 3 binaries having 
orbital parameters such that the light deflection amounts to be at least 1 /^as; see FIG. HI 

In Section IVII a finite time of observation of 5 years (Gaia mission time) has been con- 
sidered, which considerably reduces the total number of relevant binaries. Clearly, this case 
is of more practical importance, since a restricted time window of observation is in more 
agreement with reality than the first scenario. By taking into account the probability to 
find the system in the ideal astrometric position E = where the light deflection becomes 
maximal, it has been found by evaluating the corresponding integral (|27|) that there is not 
any of the relevant binary systems in the ideal position E = during Gaia mission time; 
see FIG. EJ Thus, while in principle a few binaries will have a significant light deflection, 
the effect could not be detected due to the restricted finite time window of observation. 

Furthermore, the special case of resolved binaries has been considered in Section IVlIl The 
astrom etric instrumentation of the ESA cornerstone mission Gaia, see e.g. Perryman, et al. 



2001| , has been considered in some detail in order to decide whether or not this subtle effect 
of light deflection can be observed. Two conditions for resolved binaries were presented in 
Eqs. pUj) and (1411) for such special kind of binary systems. It has been shown, however, 
that even for the Gaia mission, which is an outstanding milestone of progress in astrometry, 
such binary systems must have rather extreme orbital parameters in order to reach todays 
level of detectability, i.e. on microarcsecond level. The existence of such exotic binaries is, 
however, highly improbable. 

In summary, the main results are presented by the inclination formulae ( fl3|) and its 
simplified version ( Tl9|) . the stringent conditions ( 120]) and f[2~Tj) and by the diagrams FIG. H] 
and FIG.[5j Accordingly, one comes to the conclusion that the detectability of light deflection 
in binary systems reaches the technical limit of todays high precision astrometry and might 
be detected only in case of a very few and highly exotic binary systems. It is, however, 
very unlikely that such extreme binaries might exist. It seems that the detection of the 
light deflection effect in binary systems needs an astrometric accuracy of better than about 
0.1/zas. Thus, only astrometric missions of the next generation can accept the challenge to 
detect this subtle effect of relativity. 
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Appendix A: Two-body problem 

The calculations in this Appendix follow mainly Landau & Lifshitz |l976 |. Consider two 
massive bodies, one component having a mass M A and spatial coordinate r A , and second 
component with a mass Mb and spatial coordinate t\b, respectively. They orbit around 
their common center of mass Tcms, 

r CMS = „, ] (M A r A + M B r B ) . (Al) 
Ma + Mb 

The Lagrangian C of the two-body problem is given by 

C-^rl + ^rl-U(\r A -r B \), (A2) 

where U is the potential. With the aid of relative coordinate r ab = r 'a — r B and reduced 
mass M = Ma Mb I [Ma + Mb), the two-body problem can be transformed into into one- 
body problem, 

C = y r\ B - U(t A b) ■ (A3) 

Polar coordinates (r A B, <fi) yield 

C = I (M t\ b + r\ B 2 ) - U (tab) ■ (A4) 

The orbital angular momentum L is conserved 

L = M r AB = const , (A5) 
by means of which one obtains for the total energy of the two-body system the expression 

~M 1? 

E = T^ + IW7T B + U(rABh (A6) 



From Eq. (1A6|) one deduces 



*«,= (^lE-U(r AB )]-=^-J , (A7) 



and from Eq. flA7[) one obtains up to an integration constant 

* = / dr AB {^[E- U (^ab)} ~ =^v- ) . (A8) 



L 2 ^- 1/2 



M Z r\ B , 



-1/2 



dr AB ^- [2M[E-U (tab)} - ^~ ) , (A9) 

r AB \ r AB J 



where in the second relation Eq. flA5[) has been used; note that flA9[) is the relation between 
tab and (ft and is called orbital equation. The Eqs. flA8[) and ( 1A9j) are the general integral 
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solutions of a two-body problem. In order to integrate these Eqs. (1A8j) and (IA90 one has to 
specify the potential U. In case of Kepler problem one has 

U(r) = with a = GM A M B . (A10) 

tab 



The Eq. (IA9j) can be integrated and yields 

^ = aICCOS (^ L a j(^2M g+ T L / » j , (All) 

where the axes are chosen such that the mentioned integration constant vanishes. Further- 
more, by introducing the eccentricity e (possible values of eccentricity are between < e < 1; 
e = corresponds to a circular orbit), 

e = V + ? M\ M% ) • (A12 > 
the solution (IA110 can be written as 

1 L 2 

= = l+ecos0. (A13) 

r AB 1 MM A M B 

Note the expressions of semi-major axis A and semi- minor axis B, 

A = —SL , (A14) 

(1 -e 2 ) 7MM A M B ' V ' 

B = L l_ . (A15) 

VT^-e? 7 M M A M B v 1 

To solve the integral (IA8j) . one substitutes 

r^s - A = -A e cos E , (A16) 

where E is called eccentric anomaly. Then, one obtains for the integral in Eq. HA8|) the 
expression 

t= Um^mJ ' 72 /^' 1 — £ >- < A17 > 

and the solution is given by 

/ ,4 3 \ 1/2 
t=[—— — - (E-esinE), (A18) 

where the integration constant vanishes, i.e. the particle at t = is in periastron. The 
Eqs. (IA13I) and (|A18[) are the general solutions of two-body problem. They can be rewritten 

as 

r AB = A{l-e cosE) , (A19) 

(E-e sinE) . (A20) 



A 3 ^ 1/2 



v7 [M A + M B ) 
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In case of ellipse, E = in periastron, E = ir in apastron, and for a complete orbit E runs 
from E = 0toE = 27c. Thus, one obtains for the orbital period the expression 



T = 2vr 



A 3 



1/2 



j(M A + M B )J 

Note, the solution r in Cartesian coordinates, x = tab cos and y = tab sin 0: 



(A21) 




y 



1/2 



sin E . 



The coordinates of the bodies A and B, i.e. their orbits, are given by 

tab 



r A = r CMS + 



1 + 



Ma 



Tb = r C MS 



tab 



1 + 



Ma' 

M A 



(A22) 
(A23) 
(A24) 

(A25) 
(A26) 



Accordingly, the geometry of the orbit is determined by two orbital parameters: semi- major 
axis A and eccentricity e. In order to know the position of one celestial body, either compo- 
nent A or component B, two additional orbital parameters are needed, namely orbital period 
T and true anomaly v. A geometrical representation of the coordinates of the components 
of a binary star is given in FIG. [7] for the case of M A = 1.5 M , Mb = 1.0 M & , e = 0.5, A = 
2 AU. 
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FIG. 7: Geometrical representation of the coordinates of a binary star. In the example considered, 
the masses are Ma = 1.5 M & and Mb = 1.0 M & , respectively. The semi-major axis of the binary 
system is chosen A = 2 AU and eccentricity is taken e = 0.5. The coordinates of mass center are 
r CMS = 0- The massive bodies A and B are always in opposition to each other. 
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Appendix B: Derivation of Eq. (I10|) 



For the inclination formula the impact of eccentricity on light deflection is neglected, thus 
e = 0, implying that u> = is taken. Then, for the vectors from massive body to observer 
X\ and from massive body to source Xq, one has 

/ sin i — ei cos E \ 



Xi = r 



—ex smE 



(Bl) 



COS I 



/ cosE\ 



x = -A 



where the small parameter 



sinTJ 

V o J 

A TUB 



(B2) 



< 1 



r mA + rriB 

has been introduced. From Eqs. (1B1I) and (1B2I) one obtains for vector k 
R = xi — x , the expression 

( sin i + e 2 cos E \ 

1 



1 + 2 e2 sin i cos -E + e 2 



62 sinE 



cost 



where the small parameter 



A m,A 

6 2 = < 1 



r mA + rriB 

has been introduced. Note that (1B2j) and (1B4j) yields 

d = |fe x x \ = A |cos i| (l + O (e 2 )) . 
Using (1B1I) - (1B5j) . the generalized lens equation d^J) reads 



1 1 



2 Ti T 2 



V 2 + 8 (T + Ti - ei ) T 2 - 



.4 



where Wi = ^yl - T 2 and 



T n = sin z cos E . 



Ti 



1 — 2 ei sin z cos £? + ef 



1 + 2 e 2 sin i cos TJ + e 2 



(B3) 
R/R, where 

(B4) 



2 • 



(B5) 
(B6) 

(B7) 

(B8) 
(B9) 
(BIO) 
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By series expansion one obtains up to terms of order O 



A mA 

f \ T T 



r = \ [ ^2+8^(1+^)-^ | . (Bll) 



.4 



where W% = — Vl — w 2 with w = sin 7 cosE. The minimal and maximal light deflection 

r 

angle are 

¥?min - ¥>(l = |'^ = 7r ) =°. ( B12 ) 



^max = <^^ = -,^ = 0j =2^-— . (B13) 

In this study, the maximal possible light deflection effect is of interest. Accordingly, two 
configurations are relevant: 



1 l A 2 m A A 
m(E = 0) = - \ — cos 2 7 + 8 (1 + sinz) Icoszl I , (B14) 

2 \ V f t r r 



up to order O J™ which is just Eq. (ITU]) , and 



7T 



1 I A 2 m A A 

<p(i = ±)=-[ \j— sin 2 £ + 8— -(l + cos£)--|sin£| ) , (B15) 



(A / 777. A I 
— \ . Furthermore, it is useful to take into account only astrometric 
r V t r I 

7T 

positions with < E < — , because otherwise the light deflection is for sure negligible. 

Appendix C: Derivation of Eq. (I13|) 

From (lB14j) one obtains 

A , A 2 A 2 9 m A , i , 

2yj + — cos 7 =— cos 2 7 + 8 1 + sihj . (CI) 

ty J ty Zi tyt ty 

From (1C1D one obtains 
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Eq. (IC20 represents an quadratic equation for the expression | sinz |, which has the following 
both solutions for the inclination i: 




sinz = — ~ ±\ -. — q\ , (C3) 



where 

8m 2 A — 4mr 2 ip 2 
P = A{r 2 p 2 + 4m 2 ) ' ( ' 

A 2 r 2 p 2 + Am Ar 2 p 2 — Am 2 A 2 — r 4 p 4 
Q = A 2 (rV 2 +4m 2 ) ' ( ^ 

Eq. f lC3j) represents two solutions, however only the one with the plus-sign is valid. This can 
be shown as follows. For the value % = ^ the light deflection has to be p = <p m3X = 2y/mA/r, 
according to (|T2l) . By inserting (/3 max in Eqs. (104[) and (I05[) one obtains p = —2 m/ (A + m) 
and q = —{A — m) / {A + m). If one inserts i = S for p and g into Eq. (1031) one obtains the 
relation 



in- ^ A — m m ^ A (C6) 



A + m ^(A + m) 2 A + m A + m A + m 



Obviously, relation ( 1C6|) is only correct for the upper sign. A very similar proof can also 
be done using <£> min which also yields that the upper sign is the correct and only solution. 
Therefore, the inclination formula is given by (note, that in the region under consideration 
sini = sin (jr — i)) 



. p p 2 \ 7T 

arcsin | — - + — <?l f° r < i < — , 

(C7) 

p \ 7T . 

7T — arcsin | — 2 V ^ — / 2 

For the complete region < i < ix one obtains for the inclination formula the following 
expression: 



7T 




i 


= arccos 


2 





p p 2 



2 V 4 g| ' (C8) 



where p and q are given by Eqs. (1C4j) and ( 1C5[) . 



Appendix D: Derivation of Eq. ([3]) 



In this appendix, some basic steps of the calculations of Klioner, et al. |2003 ] are pre 



sented. A scheme of the light propagation in a binary system is presented by Fig. |5J The 
light signal from component B, considered to be the light source, is deflected by component 
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A, considered to be the massive body which curves the space-time. The vector x\ points 
from the mass center of massive body to the observer, and vector Xq points from the mass 
center of massive body to the source; R = x% — x and its absolute value R = \R\, and 

m = — — is the Schwarzschild radius of massive body; the explicit label A is omitted. 

c 2 

Furthermore we define the impact vector d = kx (x\ x k) is defined, and its absolute value 
d=\d\, see also FIG.EJ 



According to Klioner & Zschocke [2010[, the transformation of k to the unit tangent 



vector n of light-trajectory at observer is in standard post-Newtonian approach given by 

k X (Xn X £Ci) _ / 9 \ _ , 

n = k-2m — K — ^— + O (m 2 ) . (Dl) 

X\ [XqXx + Xq ■ Xi) V ' 

Note, the PPN parameter of parameterized post-Newtonian formalism, which characterizes 
a possible deviation of the physical reality from general theory of gravity, is set equal to 1 
here for simplicity. This expression is valid as long as d ^> m, but diverges for d — >■ 0. Thus, 
it is not valid for all possible binary configurations, instead one has take care to consider 
only such astrometric configurations with d ^> m. By means of fIDlj) one obtains for the 
light deflection angle (p, i.e. for the angle between n and k, the expression 

TTl IP 

cp = 2 — tan - , (D2) 

where ^ _|* ^os = ^ an z ' Xl = r ^ ^ as ^ een use< i, and ip is the angle between x 
and x±. The expression (1D2I) diverges for if) — > tt, which corresponds with the mentioned 
divergence of (1D1I) for d — > 0. Obviously, if) < i + S (from Eq. ( |28i) it is clear that eccentric 
anomaly E of binary system should be very close to zero for the light deflection effect to be 
observable at the level of microarcsecond, i.e. one actually could even assume if) ~ i + S) 
and one obtains 

m (i 7r\ m /7r i\ ,^ , 

where tan (a + IjJ = — cot a has been used, cot a = tan" 1 a, and the asymmetry of function 
cot a. From ( 1D3|) one obtains 



/ M uas pc\ ,„ , 

< 2 arctan 0.0197 — — — , D4 

KMS V M & V 9 r J 



m M 

where = — — has been used (recall M is the mass of component A and M is the 

m Mq 

solar mass), and the numerical values m ~ 1.476 x 10 3 m, 1 /ias ~ 4.848 x 10~ 12 rad and 

1 pc ~ 3.086 x 10 16 m have been inserted, so that — ~ 0.0197. The simplified inclination 

/xaspc 

formula flD4j) has first been obtained by Klioner, et al. [20031 ]: here the fact is noticed, that 
due to the divergence of the post-Newtonian solution ( ID II) for d — Y 0, which corresponds to 
ib — > 7r, the applicability of ( ID4j) is restricted by the condition d 3> m. Using d = A |cosi| 
one obtains the validity condition for the applicability of (ID4j) : 

% ^> arcsm — . (Do) 

2 kms A 
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Equation f]D4|) provides a relation between inclination % and light deflection ip for a binary 
system characerized by it's distance r from the observer and it's stellar mass M; note that 
(1D4p agrees with the simplified inclination formula given by Equation ( 1191) . 

Appendix E: Probability Distribution 

Assume a probability distribution of a quantity x is given by f(x). The probability P, 
to find a system in the interval Xi < x < Xj + Ax is given by 

/ cfe/(*) 

P(x, <x<x, + Ax) = , (El) 

dz f(z) 



where the region of validity of probability distribution /(x) is given by x m i n and x max . In 
the infinitesimal limit Ax — > dx, one obtains by series expansion the following explicit form 
for the here used probability distributions: for a power law f(x) ~ x~~ a with a^l one finds 



and for a logarithmic law f(x) ~ x 1 one has 



1 / T ^ _1 



/(x) = - ln= . (E3) 

X \ ^min / 

The normalization is J f(x) dx = 1 and the averaged value x = / f(x) xdx. 



